Measuring Gaussian Quantum Information and Correlations Using the Renyi Entropy of Order 2 
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We demonstrate that the Renyi-2 entropy provides a natural measure of infoimation for any multimode Gaus- 
sian state of quantum harmonic systems, operationally linked to the phase-space Shannon sampling entropy of 
the Wigner distribution of the state. We prove that, in the Gaussian scenario, such an entropy satisfies the strong 
subadditivity inequality, a key requirement for quantum information theory. This allows us to define and analyze 
measures of Gaussian entanglement and more general quantum correlations based on such an entropy, which 
are shown to satisfy relevant properties such as monogamy. 
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In quantum information theory, the degree of information 
contained in a quantum state p is conventionally quantified 
via the von Neumann entropy of the state S(p) = -tr(plnp), 
that is the direct counterpart to the Shannon entropy in classi- 
cal information theory [Ij]. The most fundamental mathemati- 
cal implications and physical insights in quantum information 
theory, ranging from the Holevo bound |2i| all the way to the 
whole genealogy of quantum communication protocols ylQ], 
rely on a key property satisfied by th_e_von Neumann entropy, 
the strong subadditivity inequality | 



S{pab) + SipBc) < S(pabc) + S(pb) , 



(1) 



for an arbitrary tripartite state pabc- The strong subadditivity 
inequality implies in particular that the mutual information 



I(Pa:b) = S(pa) + S(pb) - S(pab) , 



(2) 



which measures the total correlations between the subsystems 
A and B in the bipartite state pab, is always nonnegative. How- 
ever, in classical as well as quantum information theory, sev- 
eral other entropic quantities have been introduced and stud- 
ied. In particular, Renyi-a entropies ^ are an interesting 
family of additive entropies, whose interpretation is related 
to derivatives of the free energy with respect to temperature 
1^, and which have found applications especially in the study 
of channel capacities [^{|, work value of information lUOll . and 
entanglement spectra in many -body systems jllll . They are 
defined as 



>S,(p) = (!-«)-' In tr(p"). 



(3) 



and reduce to the von Neumann entropy in the limit ff — > 1 . 
Renyi entropies are powerful quantities for studying quantum 
correlations in multipartite states. For a bipartite pure state 
Pab, any of the entropies in Eq. Q evaluated on the reduced 
density matrix of one subsystem only, say pA, is an entan- 
glement monotone uM, dubbed Renyi-a entanglement lOl . 
Any such measure can be extended to mixed states via con- 
ventional convex roof techniques UM . Using the a = 2 in- 
stance, S2(p) = -lntr(p-), the ensuing Renyi-2 measure of 
entanglement has been defined and proven to satisfy the im- 
portant 'monogamy' inequality lITsl II6II for multiqubit states 
ilSlUTIl . Nonetheless, Renyi-a entropies for a ?t 1 are not in 



general subadditive ^; this entails, e.g., that if one replaces 
S by Sa^i in (|2|, the corresponding quantity can become neg- 
ative, i.e., meaningless as a correlation measure; see lUSll for 
explicit two-qubit instances of this fact when a -2. 

In this Letter, we focus our attention on multimode quan- 
tum harmonic oscillators, and analyze the informational prop- 
erties of general Gaussian states as measured by the Renyi- 
2 entropy. Gaussian states constitute versatile resources for 
quantum communication protocols with continuous variables 
lUOluOll and are important testbeds for investigating the struc- 
ture of quantum coiTelations 1121 1 , whose role is crucial in 
fields as diverse as quantum field theory, solid state physics, 
quantum stochastic processes, and open system dynamics. 
Gaussian states naturally occur as ground or thermal equi- 
librium states of any physical quantum system in the 'small- 
oscillations' limit i22t i23i1 , and can be very efficiently engi- 
neered, controlled and detected in various experimental se- 
tups, including light, atomic ensembles, trapped ions, nano- 
/opto-mechanical resonators, and hybrid interfaces thereof 
II20I1 . Given the special role played by Gaussian states and 
operations, which are formally 'extremal' with respect to sev- 
eral quantum primitives ll23ll . a Gaussian-only theory of quan- 
tum information is actively pursued |[24ll25ll . In particular, a 
key conjecture which — despite its disproof for general chan- 
nels 1I26I1 — still stands within the Gaussian scenario, is the so- 
called minimum output entropy conjecture for bosonic chan- 
nels ll24ll25ll27ii28il . Interestingly, the conjecture is verified 
(for all phase-insensitive channels) if using Renyi-a entropies 
with a > 2 1I29I1 . but to date is resisting an analytical proof for 
a ^> \ II30I1 . This somehow raises the question whether the 
von Neumann entropy is indeed the most natural one within 
the Gaussian scenario Olll . 

Here we show that the Renyi-2 entropy should be regarded 
as a specially meaningful choice to develop a Gaussian theory 
of quantum information and correlations. We prove that the 
Renyi-2 entropy S2 satisfies the strong subadditivity inequal- 
ity ([Til for all Gaussian states, is operationally linked to the 
Shannon entropy of Gaussian Wigner distiibutions (providing 
a natural measure of state distinguishability in phase space), 
and can be employed to define valid measures of Gaussian en- 
tanglement JSZ, i33J and general discord-like quantum corre- 
lations IBJI - OSll . whose monogamy properties lITsl [la, I39l - l42ll 
we investigate in detail. In particular, we obtain a truly bona 



fide measure of genuine tripartite entanglement for three- 
mode Gaussian states, based on 82- Our study allows us to 
explore within a unified framework the various facets of non- 
classicality in the Gaussian realm. 

Renyi-2 entropy for Gaussian states. — We consider a n- 
mode continuous variable system; we collect the quadrature 
operators in the vector R - {q\,p\,q2, Pi, ■ ■ ■, Qm PnV 6 R-^" 
and write the canonical commutation relations compactly as 
[Rj,Rk] - i{of'")ji^ with o) - ( j A being the symplectic ma- 
trix. A Gaussian state p is described by a positive, Gaussian- 
shaped Wigner distribution in phase space. 



Wp(^) 



1 



n" ^/deVy 



exp(-^'7"'^) 



(4) 



where ^ e R^", and 7 is the covariance matrix (CM) of ele- 
ments yj^it - ii[p{Rj,Rt]+], which (up to local displacements) 
completely characterizes the Gaussian state p. 

Let us now evaluate different entropic quantities on the 
Gaussian state p. All the measures defined below are invari- 
ant under local unitaries, so we will assume our states to have 
zero first moments, (R) - 0, without loss of generality. The 
purity is easily computed as trp^ = {2n)" f ,„ W^i^) d^"^ - 

(dety) 2 . Hence the Renyi-2 entropy of an arbitrary n-mode 
Gaussian state is 



Slip) = I Mdstj)., 



(5) 



ranging from on pure states (det/ - 1) and growing un- 
boundedly with increasing mixedness (i.e., temperature) of 
the state. The von Neumann entropy is instead a more com- 
plicated function that depends on the local temperatures of 
all the « normal modes of p, that is, on the full symplec- 
tic spectrum of 7 1I21I l3ll l43l 14411 . On the other hand, since 
the Wigner function is a valid probability distribution provid- 
ing a fully equivalent description of p, we can alternatively 
compute the Shannon entropy of Wp to obtain an alternative 
quantifier of the informational content of any Gaussian state. 
Such an entropy has a clear interpretation in terms of phase- 
space sampling via homodyne detections ||45I1 . The contin- 
uous (Boltzmann-Gibbs-)Shannon entropy for a probability 
distribution P{x) is defined as H{P) = - f Pjx) In P(x) dx 
Il46i1 . It can be shown that II47I1 (see the Appendix 14811 for 
a derivation) 

H(Wp) = -JjVp(^)MWpmd^"^ 

= Slip) + nil + Inn). (6) 

Interestingly, we see that this sampling entropy coincides 
(modulo the additional constant) with the Renyi-2 entropy 
rather than the von Neumann one for Gaussian states. We can 
now take an extra step and introduce a 'relative sampling en- 
tropy', to quantify the phase-space distinguishability of two 
«-mode Gaussian quantum states pi and p2 (with CMs 7; 
and 72, respectively), defined as the relative Shannon entropy 
(also known as KuUback-Leibler divergence j^ |a]) between 



their respective Wigner distributions Wp, and Wp.., yielding 



HiWpMp^) ± f Wp,(^ln( 



Wp,i^ 



- H'"(S^- 



WpM) 



d'"€ 



(7) 



■ « . 



Let us evaluate Eq. (|7]i when pi = pab is a generic Gaussian 
state of a composite system, partitioned into two subsystems 
A and B (of n^ and n^ modes respectively, with «^ + n^ - «), 
andp2 =Pa®Pb is the tensor product of the two marginals of 
Pab- Writing the CMs in block form. 



7\ = Jab 



7a 



^ab 



^AB 



7b 



7i=7a®7b, 



(8) 



we have tr(7i72M — 2inA + hb) - 2«, which entails 

HiWpJ\Wp,^p,) = HiWp,) + HiWp,) - HiWp,J 

1 ^^/ det7^det7jj 

2 I del Jab 

= SlipA) + SlipB) - SlipAB) 
- IliPA-.s)- 



(9) 



The above equation defines the 'Gaussian Renyi-2 (GR2) 
mutual information' I2 for an arbitrary bipartite Gaussian 
state Pab- It follows that l2iPA:B) ^ 0, as it coincides ex- 
actly with the Shannon continuous mutual information of the 
Wigner function of pab, which is positive semidefinite. The 
expression in Eq. ^, which is analogous to the more fa- 
miliar von Neumann one [Eq. Q], has thus a precise oper- 
ational interpretation as the amount of extra discrete informa- 
tion (measured in natural bits) that needs to be transmitted 
over a continuous variable channel to reconstruct the com- 
plete joint Wigner function of pab rather than just the two 
marginal Wigner functions of the subsystems 11451 14611 : in 
short, l2iPA:B) measures the total quadrature correlations of 

PAB- 

We can then enquire whether a more general property such 
as the strong subadditivity inequality ([TJ — which would in- 
dependently imply the nonnegativity of J2 — holds in general 
for S2 in the Gaussian scenario. Let Pabc be ^ tripartite Gaus- 
sian state whose subsystems encompass arbitrary number of 
modes. Writing its CM in block form as in Eq. (|8), and using 
the definition (|5j, we have the following 

Theorem 1 The Renyi-2 entropy S2 satisfies the strong sub- 
additivity inequality for all Gaussian states Pabc> 



S2ipAB) + SlipBc) - SlipABc) " SlipB) 

^iJdet7ABdetrBc\^ 
2 \det7^Bcdet7B/ 



(10) 



Proof. The result follows by applying a particular norm 
compression inequality to the CM Yabc- Given a positive 
Hermitian matrix A e M„,, and given any two index sets 
a,f5 Q N - {!,..., m), the Hadamard-Fisher inequality ll49ll 



states that det Aq-u/j det Ao-n^s ^ detA^, det A^. Recal ling that 
any CM y^^^ i^ ^ positive real symmetric matrix IdOII . the 
claim follows upon identifying a with the indices of modes 
AB and p with the indices of modes BC. ■ 

Beyond its apparent simplicity. Theorem [T] has profound 
consequences. It yields that the core of quantum information 
theory can be consistently reformulated, within the Gaussian 
scenario 11241 12511 . using the simpler and physically natural 
Renyi-2 entropy in alternative to the von Neumann one. In the 
rest of this Letter we will focus on defining GR2 quantifiers 
of entanglement and other correlations for Gaussian states. 
Gaussian Renyi-2 measures of correlations. — The GR2 en- 
tanglement £2 can be defined as follows. Given a «-mode bi- 
partite Gaussian state pab with CM 7^^, 



&2(Pa:b) - 



inf 

"'AB- 

^^AB^yAB- 
dcltr^B = l 



w®" 



- In (det a- a) ■ 



(11) 



For a pure Gaussian state pab = I'Aab)(i/'abL the minimum 
is saturated by ctab = 7aB' so that &2{4'A:b) - -SiifA) - 
^ln(det7^), where 7^ is the reduced CM of subsystem A. 
For a generally mixed state, Eq. (fTTT i — where the minimiza- 
tion is over pure n-mode Gaussian states with CM (Tab smaller 
than y^g — amounts to taking the Gaussian convex roof of the 
pure-state Renyi-2 entropy of entanglement, according to the 
formalism of lull 13211 . Closed formulae for S2 can be obtained 
for special classes of two-mode Gaussian states |i33i i51i |52i] 
exploiting the same procedure adopted for the Gaussian en- 
tanglement of formation II32II48II . Like the latter measure (and 
all Gaussian convex-roof entanglement measures), it follows 
from the results of ll32ll that the GR2 entanglement is in gen- 
eral monotonically nonincreasing under Gaussian local oper- 
ations and classical communication, and is additive for two- 
mode symmetric Gaussian states. 

We can also introduce a 'GR2 measure of one-way classi- 
cal correlations' in the spirit of Henderson and Vedral IBSll . 
We define J'2{pa\b) as the maximum decrease in the Renyi-2 
entropy of subsystem A, given a Gaussian measurement has 
been performed on subsystem B, where the maximization is 
over all Gaussian measurements — i.e. those that map Gaus- 
sian states into Gaussian states OTllSSll . Any such measure- 
ment on, say, the «B-mode subsystem B - (Bi . . . B„g), is de- 
scribed by a positive operator valued measure (POVM) ll53ll 
of the form n^d/) = n-"nU%i WB,(r]j)]A^[n%i wlir^j)] 
where Wsirij) - exp(rijb . — ri*bj) is the Weyl operator, bj - 

(qBj + ipB,)/ V2, TT'"" jnB{n)d^"''il = 1, and A" is the density 
matrix of a (generally mixed) n^-mode Gaussian state with 
CM r^ which denotes the seed of the measurement. The con- 
ditional state pa\jj of subsystem A after the measurement TlBitl) 
has been performed on B has a CM 7J independent of the out- 
come I] and given by the Schur complement |54[| 



fA^7A-^AB(7B+T'^r'?lB^ 



(12) 



where the original bipartite CM Yab of 'he «-mode state pab 
has been written in block form as in dHJ. We have then 



J2(Pa\b) 



sup - In 



det 7^ 
detT-n 



(13) 



The one-way classical correlations J'2{pb\a), with measure- 
ments on A, can be defined accordingly by swapping A «-» B. 
We can now define a Gaussian measure of quantumness of 
correlations based on Renyi-2 entropy. Following the seminal 
study by Ollivier and Zurek ll34ll . and the recent analysis of 
Gaussian quantum discord using von Neumann entropy 11371 
13811 . we define the 'GR2 discord' as the difference between 
mutual information (|9j and classical correlations (fTJt . 



D2(PA\b) - I2(PA:b) - J'2(PA\b) 

1 Cdet7Bdet7n^ 
= mf - In 

rn 2 l^ dety^B 



(14) 



For the case of A and B being single modes, that is, pab being 
a general two-mode Gaussian state, closed formulae can be 
obtained for Eqs. (113114b thanks to the results of 1371. We re- 
port them in the Appendix [48] for completeness. We remark 
that Theorem [1] is crucial to guarantee the nonnegativity and 
the faithfulness of the GR2 discord 1i55h57i1 . Notice also that 
\l2ipA:B) = J2(Pa\b) = •2)2(P/1|b) = ^2(Pa) for purc bipartite 
Gaussian states pab- 

A trade-off relation between the entanglement £2 and the 
classical correlations J2 can be written for arbitrary tripartite 
pure Gaussian states pabc 1-371], following Koashi and Win- 
ter |59]. One can essentially exploit the fact that all possible 
Gaussian POVMs on B induce all possible Gaussian pure- 
state decompositions of the subsystem AC (see also 11611 ). 
which implies 



S2(Pa) = J2(PA\b) + £>2(PA:c) ■ 



(15) 



This relation can be manipulated to express the "conservation" 
of different types of correlations in a generic pure tripartite 
Gaussian state pabc^ along the lines of 1I6OI1 . 

We now look at monogamy properties of the GR2 mea- 
sures. For an entanglement monotone E and a «-partite 
state PA|/1i...a„, the monogamy relation (choosing party Ai 
as the focus), which constrains the distribution of bipartite 
entanglement among different splits, can be written as lITsIl 
E(pAi:A2...AArzX]=2 E(pa,:A,) > . The Renyi-2 entanglement 
measure Ill3lll7ll . as well as the tangle (squared concurrence) 
ITsl Uql . satisfy this inequality for general «-qubit states. A 
Gaussian version of the tangle (based on squared negativ- 
ity) has been defined that obeys the inequality for all «-mode 
Gaussian states ll40ll. We now show that £t does too. 



Tlieorem 2 The GR2 entanglement defined in Eq. ( 1771 ) is 
monogamous for all n-mode Gaussian states pa,a2...a„, 



£2(PA,:A2..a) - Ii"=2 £'2(PAi:Aj) > , 

where each Aj comprises one mode only. 



(16) 



Proof. The structure of the proof follows closely the one for 
the tangle of n-qubit systems [1(S]. It suffices to prove the 
inequality for tripartite Gaussian states PA1A2A' where A3 com- 
prises « - 2 modes (with n arbitrary), as iterative applications 
to fi2(PA,:A') would then imply Eq. ( fT6l ). It is further enough 
to prove the inequality on pure states, as it would then ex- 
tend to mixed ones by convexity lUSl [iq, I39ll40tl . Exploiting 



the phase-space Schmidt decomposition 112 IL I43L 16111 . when 
PA1A2A' is pure, the state of subsystem A'^ is locally equiva- 
lent to a 2-mode state, tensored by « - 4 irrelevant vacuum 
modes. The problem reduces to proving that &2{Pa:Bc) ^ 
&2{Pa:b) + S'liPA-.c) for an arbitrary pure Gaussian state pabc, 
where A and B are single modes, while C groups two modes. 
Noting that SoiPA-.Bc) = 'S2(pa), and exploiting Eq. (flSl l, we 
see that Eq. (fTSl l is verified if one establishes that 



J2{PA\b) > SliPA-.B) 



(17) 



holds for all (mixed) two-mode Gaussian states p^^. The in- 
equality ([TtT i is proven in the Appendix |48], which concludes 
the proof of the Theorem. ■ 

Let us analyze in detail the case of Pa^a^a^ being a pure 
3-mode Gaussian state, whose CM is characterized up to lo- 
cal unitaries by three parameters (local symplectic invariants) 
Uj > 1, with SiipAj) = Inoj (j = 1,2,3) HH. We define 
the residual entanglement emerging from the monogamy in- 
equality as £;2(pa,:A,:A3) = £'2(PA.:A-.aJ-£'2(PA.:A.)-£'2(PA,:A,)- 

This quantity, which can be calculated exactly 11391 [5211 . de- 
pends in general on the focus mode (say Ai) chosen for the 
decomposition of the bipartite entanglements. Remarkably, 
we find that for all pure three-mode Gaussian states which are 
fully inseparable and display entanglement in all global and 
reduced bipartitions (£2(Pa,:A,) > 0, &2(PAr.AjAt) > 0, V/ ?i 
j + k), the residual GR2 entanglement S2{pA,:A^:Ai) is invari- 
ant under mode permutations, thus representing (to the best 
of our knowledge) the first and only known intrinsically bona 
fide measure of genuine tripartite Gaussian entanglement ll62ll . 
We report its explicit formula here, while a derivation is pro- 
vided in H, 



£2(Pa,:A,:A3) = In(8aifl2a3) (18) 

with 6 = n^^ v=o ((ai + (-l)^fl2 + {-iya3? - 1). The formula 
( fTSl l holds when la, - a,| + 1 < oj. < (a^ -H a^ - 1)2. 

We finally remark that, importantly, the GR2 discord 
Eq. (fl4] i also turns out to be a monogamous measure of 
quantum correlations for arbitrary pure three-mode Gaussian 
states p J, /I ,,4, (unlike the von Neumann entropy-based discord 
illlllEll])- By using Eq. O, one finds (see also iHIMl) 

D2(PA^:A2:A,) - S^liPA.JA.A,) - i)2{PAMl) ' ^liPAMi) = 

&2{pAi\A^:Ay)- In Other words, the residual tripartite discord 
equates the residual tripartite entanglement for pure pa^AiA^- 
This extends to the multipartite case the equivalence between 
entanglement a nd g eneral quantum correlations valid for pure 
bipartite states Il34ll64ll65ll . and places D2 as the only known 



measure of quantumness beyond entanglement in continuous 
variable systems that fulfills monogamy |;66;1 . 

Conclusions. — In this Letter we planted the seeds for a full 
Gaussian quantum information theory ll25ll using the Renyi-2 
entropy 82- This is possible thanks to the fact, proven in The- 
orem[ri that such an entropy satisfies the strong subadditivity 
inequality for arbitrary Gaussian states of quantum harmonic 
systems. We employed S2 to define valid measures of entan- 
glement, total, classical, and quantum correlations, highlight- 
ing their properties. The Renyi-2 mutual information is inti- 
mately related to Wigner distribution sampling by homodyne 
detections in phase space. The residual Renyi-2 entanglement 
measure allows for a quantification of genuine tripartite entan- 
glement in three-mode fully inseparable pure Gaussian states, 
which is invariant under mode permutations. We argue that 
the measures defined in this Letter should be adopted as priv- 
ileged tools to address the quantification of relevant correla- 
tions in Gaussian states. A very recent application to relativis- 
tic quantum information has been reported ll68ll . 

By conception, this work has been biased towards Gaus- 
sian states and operations. However, the approach pursued 
here can be extended to arbitrary, even non-Gaussian «-mode 
states p of continuous variable systems, provided one chooses 
the Wehrl entropy H to quantify the informational content 
of p. Such an entropy is operationally associated to phase- 
space sampling via heterodyne detections ll45ll as it corre- 
sponds to the continuous Boltzmann-Gibbs-Shannon entropy 
of the Husimi Q distribution of p, Q{a) - n^"{a\p\a), which is 
a valid (nonnegative) probability distribution for all quantum 
states llll. One can then adopt a distance measure between 
any two p\ and p2 in terms of the relative Shannon entropy 
between their respective Q distributions, define ensuing cor- 
relation measures, and so on. Such a formalism could also ac- 
commodate a measure of non-Gaussianity of quantum states 
1691], and one might then naturally compare Gaussian with 
non-Gaussian operations for the realization of specific tasks 
II27I1 such as extracting classical correlations IITOl ItiIi . maxi- 
mizing information-disturbance trade-off" ll72ll and performing 
optimal cloning of coherent states ||73I1 . This will be the sub- 
ject of further study. 
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Appendix A: Derivation of formulae © and Q 

For the sake of completeness, here we provide the reader with an explicit derivation of Eqs. (|6]l and O, on which the opera- 
tional interpretation of the Renyi-2 entropy is ultimately based. 
Direct substitution of Eq. (HJi into (|6]l yields 



H{Wp) 



r -^ 



Vdety 



exp(-^'r"'^) 



^'7 '^ + «ln;r+-ln(detr) 



d'"^ . 



(Al) 



The first term in the square bracket is conveniently handled by performing the integration in phase space coordinates ^ that 
diagonalize the symmetric, positive definite matrix 7 (whose eigenvalues will be denoted by jj), and by noting that, for any 
7 > 0, one has 



-f 

J^ J-o. 



exp(-7 \x^)x^dx = ^ 



(A2) 



Hence, by applying ( IA2I ) as well as the normalization of the Wigner function to ( lAU . one obtains Eq. ©I 



H{Wp) = y — + nXmi + Snip) = « + n\nn + S2ip) 



Next, the substitution of (|4]i into (|7]i leads to: 

H{Wp,\\Wp,)^-H{Wp,)+ r ^ exp(-^V^) 

Jr-" n" ydet7i 



^^y^^^ + n\nn+-\n{A&iy2) 



j2n 



d^"^. 



(A3) 



(A4) 



Once again, it is expedient to carry out the integration in phase space coordinates where y^ is diagonal (with eigenvalues 72^-, 
while the entries of y^ will be denoted by yijk)- Then, Eq. (IA2b and the normalization of Wp^ lead to (|7]i 



H{Wp,\\Wp,) ^ -\n 



1 /det72 



2n 



y\,jj 1 , /det72 



- « + > — ^ = - In , 

j-^ 272 J 2 \det7i 



-« + -tr(7i72'), 



2 \det7i 
where the last step follows from the invariance of the quantity tr(7i72') under changes of basis 



(A5) 



Appendix B: Explicit formulae for two-mode Gaussian correlations 



Standard form. — The CM 7^^ of any two-mode Gaussian state pab can be transformed, by means of local unitary operations, 
into a standard form of the type lulll 



Tab 



7a ^ab ] 




^Ib Tb J 





a c+ ^ 

fl c- 

c+ b 

c_ b ) 



(Bl) 



where a,b > 1 , \\cp- - l) (b^ - l) - 2c_c+ - abc^ + c-_ i-ab + c\\\ > 0, and we can set c+ > |c_| without losing any generality. 

These conditions ensure that the uncertainty relation y^s ^ '^®^ is verified, which is a bona fide requirement for the CM y^B to 
be associated with a physical Gaussian state in a two-mode infinite-dimensional Hilbert space 1I21I1 . Recall that for pure Gaussian 

states, b - a,c+ - -c_ = Va- - 1. 



All the formulae presented in the following will be written explicitly for standard form CMs for simplicity. However, they 
can be recast in a locally invariant form by expressing them in terms of the four local symplectic invariants of a generic two- 
mode Gaussian state ll44ll . /i = dety^, I2 - det/g, 73 = det^^^, U - dety^^. This is accomplished by inverting the relations 
/i - a^, h - b^, h - C+C-, I4 - (ab — c+){ab - c_) so that the {Ijf-^^ appear explicitly in the formulae below 1I21I1 . The obtained 
expressions would then be valid for two-mode CMs in any symplectic basis, beyond the standard form. 

GR2 entanglement. — For generally mixed two-mode Gaussian states p^B, the Renyi-2 entanglement measure £2(Pa:b), defined 
by Eq. (fTTl i in the main text, admits the following expression if the CM 7^^ is in standard form ll32ll33ll . 



&2{Pa:b) = :t In inf nieia, b, c+, c_) , 

2 \ee[0,2n] ' 



(B2) 



with 



tng(a, b, c+, c_) 



1 + 



c+(ab - c_) - c_ 



cos J\a - b(ab — c^) \\b — a(ab — c^)\ 



X < 



2 {ab - c^) {cp- + b^ + 2c+c_) + sin {a^ - b') 



\c+(ab — c^) + c_ 



"^ [a- b(ab - cl)] [b - a(ab - cl)] 

^ -I 

COS [label + (fl2 + ^2) c+c-_ + (( 1 - 2b^) a^ + b') c_ - ab {a- + b^--2) c+] 



^[fl - biab - c2)] [b - a(ab - cl)] 



(B3) 



The optimal minimizing Eq. ( IB3I ) can be found numerically for general two-mode Gaussian states 13211 . and analytically for 
relevant subclasses of states (including symmetric states JSlll . squeezed thermal states, and so-called GLEMS — Gaussian states 
of partial minimum uncertainty OSll '). 

GR2 classical correlations and discord. — For generally mixed two-mode Gaussian states pab, the Renyi-2 measures of one- 
way classical correlations J'2(pa\b) and quantum discord D2{pa\b), defined by Eqs. (fT3b and (fT4l) in the main text, respectively, 
admit the following expression if the CM j^b is in standard form ll37ll 

(B4) 
(B5) 



J2(Pa\b) = In a - - In inf det 7" '* , 

Z \i.(p j 

D2(pa\b) = In fc - - In ( det Jab) + :t 1" i"f '^^^ ^a " 



with A € (0, 00), (^ e [0, 2n], and 

2a\b + A){\+ bX) -aicl+ d] (2bA + A^ + l] + 2c\cIA + a{c\- c^) [a^ - l) cos(2(/j) 

det 7,'* = i^ — i^ . (B6) 

'^ 2(b^A\\+bA) 

The optimal values of A and y? minimizing Eq. ( IB6I ) can be found analytically' for all two-mode Gaussian states ll37ll . In 
particular, for standard form CMs, one gets 



inf det 7 , ' 



«(«-t) 



if [ab^cl -c\{a + bc~_)) [abV^ -c-[a + bcl)) < ; 



2\c_c+\-\/{a{h'—\)-hc-_){a(b'—l)-bcl)+(a(b^-\)-bcl)(a(b--l)-bcl)+clcl 



Otherwise. 



(B7) 



Inserting Eq. ( |B7t into Eqs. (IB4IB5I ) one gets closed formulae for the one-way GR2 classical correlations and for the GR2 
discord of general two-mode Gaussian states. 



' In general, given that <S2(Ag) is concave on tlie convex hull of the set of 
Gaussian states — having inherited the property from the Shannon entropy 
of the corresponding Wigner distributions, via Eq. (|6) — and given that 
every Gaussian state admits a convex decomposition into pure Gaussian 
states (see i28|,[37ll7Cll for more details), it follows that the optimizations in 
Eqs. U3ll4t for an arbitrary number of modes are always achieved by pure 
Gaussian seed elements, i.e., det(r2) = 1. This simplifies considerably the 
evaluation of GR2 one-way classical correlations and discord. 



Appendix C: Proof that GR2 classical correlations exceed GR2 entanglement 

Here we prove that the inequaUty 

J2{PA\b)>&i(PA:b), (CI) 

holds for all two-mode Gaussian states pab- This is a central step in the proof of the general monogamy inequality for the GR2 
entanglement measure, reported in Theorem |2]in the main text (see also HQ]). Without loss of generality, we can assume the 
CM Yj^g in standard form. 

We observe from Eqs. (IB2IB4b that setting any value of in ( IB3I ) provides an upper bound to &2, while setting any value of 
A, if in ( IB6I ) provides a lower bound to J/a, 

&2{PA:b) < -\n{mg{a,b,c+,cS)), 



JiiPMB) > lnfl-iln(detr"'-^). 



We will then set = tt (similarly to what done in 114011 ). and A - 1,^ = (corresponding to heterodyne detections on B [13811 ). 
and proceed to prove that In a - i In (det 7^ ' °) ^ 5 In {mAa, b, c+, c_)), which impUes dCll ). We then want to prove the inequality 
l+F^/G<J, where 

F = c+iab- c_) - C- + J\a- biab- c-^)\\b - a(ab - c^)\, 

label + (fl^ + b^) c+cl + ((1 - 2b^) a^ + b^) c_ - ab {a^ + b^-2) c+ 



G = 2[ab-cl)[a^ +b- + 2c^c+)- 



J = 



^j\a - b(ab — c^) \\b — a(ab - c^)\ 
a^l + bf 



(a + ab — c~)(fl + ab — c~_) 



This is equivalent to show that, defining K - K(a, b, c+, c_) = 7 - F^/G - 1, we have min(„j, f. £. ) K = Q. 

Let us search for the absolute minimum of K. This function has no singularities apart from the trivial case of two-mode 
vacuum states [a = b = l,c± - 0) for which we know that lim(nj, c+,f_)^(i,i,o,0) K = 0. Thus, we can focus on the stationary points 
of K and on the values at the boundaries of its domain. We adopt the Karush-Kuhn-Tucker method |74ll751 . Given constraints 
/ = fj{a,b,c+,C-) < and associated multipliers A, > 0, the coordinates {a*,b*,c*^,c*_,A*) associated to a minimum of K 
satisfy the following conditions: 

(a) (V, ^) L(a, b, c+, c_, A,)|(„.,iv:.c-,Ap = (v, ^) (K + 2, A,/i)|(«.,&.,<,c-,A;) = 

(b) Aifiia*, b", cl,c''_) = 0, A,- > V/ 

(c) there exists a vector x - {x,} such that (VK - 2i XiVfi)\(^a'.h'.ci,c'_) - 0. 
Conditions (a) and (b) are necessary, while condition (c) is sufficient. In this case, we have 

/i = a^ +b'^ + 2c+c_ -l-{ab- cl){ab - c'i) , 
/2 = -fl2 _ ^2 ^ 2c+c_ + l+{ab- cl){ab - cl) , 
h = -c+ - c_ . 

A bit of algebra reveals that condition (a) is verified by the solutions to the following system 

i dKdf^^dKdfi^ 

da db db da 



dK _dK ^ f(|^-|^) ^ ,^ p/i ^ df2 dfi dh 

dc- dc+ Ml " \ dc+ dc+ 5c_ dc- 

da 



(C2) 



It follows that, among the stationary points, the absolute minimum of K is again reached at a = ^ = 1 , c± = 0, yielding K - 0: it 
is immediate to verify that the conditions (b) and (c) are satisfied as well at (1, 1, 0, 0). Finally, we have to check the values of K 



at the boundaries of its domain. The relevant cases arefl-^oo,/7^oo,c+= ^ab - 1 , c+ = |c_|, c_ = 0, and somewhat tedious 
yet straightforward analysis reveals that K >Q always holds for all those cases. 

We therefore conclude that the function K is always nonnegative, which proves Eq. dClb . thereby proving that the entan- 
glement measure £2 satisfies the general monogamy inequality [Eq. (fTSI l of the main text] for all «-mode Gaussian states of 
continuous variable quantum systems. ■ 



Appendix D: Residual tripartite GR2 entanglement for pure three-mode Gaussian states 

Up to local unitaries, the CM 7^^^^^^ of any pure three-mode Gaussian state can be written in the following standard form 



7a,a,a, - 
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fl2 





< 








c- 





ai 
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a3 
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Cj 





cT 





a-i ) 



(Dl) 



where 



V[(fl,- - 1)2 - {Oj - akf-][{ai +\)~- (aj - atf] ± ^J[{ai - 1)2 - (uj + ak)^][(ai + l)^ - (aj + at)^] 



4 yjapk 



and \aj - ak\ + \ < a,- < aj + a^ — I , 



with {/, j, k] being all possible permutations of { 1, 2, 3). 

The GR2 entanglement in the two-mode reduced state with CM 7^^ is 



(D2) 



withlBJ 



8k 



1 , if at > Jfl2 -I- fl2 _ 1 ; 



8fl2' 



-I- \ 



ar — I 

k 



if aif < at < Jar + a'. - I; 
if Ok < au. 



(D3) 



Here we have set 



Ok 



\ 



2(fl2 + fl2) + (a2 - a2)2 + 1^2 _ ^2| ^(a2_^2)2 + 8(«2 + ^2) 



2(fl2 + fl2) 

P = —I + 2a[ -H 2fl2 + 2a3 -H 2fl[fl2 + 20^0^ + ^a^a-^ — a-^- a^ — a^— yd , 

6 — (-1 + fll - fl2 - fl3)(l + fll — fl2 — fl3)(— 1 + fll + fl2 — fl3)(l H- fll + fl2 - fls) 
X (-1 + fll — fl2 -(- fl3)(l + fll — fl2 + fl3)(-l -(- fll + fl2 H- fl3)(l H- fli 4- fl2 + fl3) . 



The residual GR2 entanglement, with respect to the focus mode A,, is 



S'liPAiAjAt) - ^liPArAjAt) - S'liPAiAj) - ^liPArAt) 

1 ( a- - 

2 [gkgj 



(D4) 



In general, this expression is dependent on the choice of the focus mode. However, let us consider the relevant case of a fully 
inseparable three-mode pure Gaussian state, for which entanglement is nonzero for all global splittings and for all reduced 
two-mode bipartitions, £2(Pa,:A aJ > 0, &2(pAr.A ) > 0, V{/, 7, k). In our parametrization, this occurs when ll33ll 



2i - aj\ + !<«<:< Ja^ + a^. — \ , 



(D5) 



for all mode permutations. It is immediate to see that the simultaneous verification of such a condition for all mode permutations 
imposes a^ > Uk'ik - 1, 2, 3. In this case, exploiting Eq. ( |D3t , the residual GR2 entanglement becomes 



&l(PA,:AfA,) = 2 ^" 



(6Aa]a]al^ 

I J K 



(D6) 



which is manifestly invariant under mode permutations, as reported in Eq. (fTSl l in the main text. This symmetry is broken on 
states for which some of the reduced two-mode bipartitions become separable. 



